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We calculate the masses and leptonic decay constants of the heavy vector quarko- 
nia, J/tp and T mesons at finite temperature. In particular, considering the thermal 
spectral density as well as additional operators coming up at finite temperature, the 
thermal QCD sum rules are acquired. Our numerical calculations demonstrate that 
the masses and decay constants are insensitive to the variation of temperature up 
to T = 100 MeV, however after this point, they start to fall altering the temper- 
ature. At deconfinement temperature, the decay constants attain roughly to 45% 
of their vacuum values, while the masses are diminished about 12%, and 2.5% for 
J/ijj and T states, respectively. The obtained results at zero temperature are in 
good consistency with the existing experimental data as well as predictions of the 
other nonperturbative models. Considerable decreasing in the values of the decay 
constants can be considered as a sign of the quark gluon plasma phase transition. 

PACS numbers: 11.55.Hx, 14.40.Pq, ll.10.Wx 
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I. INTRODUCTION 



Investigation of the heavy mesons canplay essential role in understanding the vacuum 
properties of the nonperturbative QCD [l|. In particular, analysis of the variation of the 
parameters of the heavy quarkonia, namely bottonium (bb) and charmonium (cc) in hadronic 
medium with respect to the temperature can give valuable information about the QCD vac- 
uum and transition to the quark gluon plasma (QGP) phase. Determination of the hadronic 
properties of the vector mesons in hot and dense QCD medium has become one of the most 
important research subject in the last twenty years both theoretically and experimentally. 
J ft/) suppression effect due to color screening can be considered as an important evidence 
for QGP 0- This suppression effect has been observed experimentally in heavy ion collision 
experiments held in super proton synchrotron (SPS) at CERN and relativistic heavy ion 
collider (RHIC) at BNL. 

Properties of the heavy mesons in vacuum have been probed widely in the literature using 
the nonperturbative approaches like QCD sum rules, nonrelativistic potential models, lattice 
theory, heavy quark effective theory and chiral perturbation theory. However, in expansion 
of most of these models to finite temperature we are face to face with some difficulties. 
QCD sum rules which is based on the operator product expansion (OPE), QCD Lagrangian 
and quark-hadron duality, is one of the most informative, applicable and predictive models 
in hadron physics [1, 0]. The thermal version of this model proposed by Bochkarev and 
Shaposhnikov has some new features at T ^ One of the new feature is the 

interaction of the particles existing in the medium with the currents which demands the 
modification of the hadronic spectral function. The other new picture of the thermal QCD 
is breakdown of the Lorentz invariance via the choice of reference frame. Due to residual 
0(3) symmetry at finite temperature, more operators with the same dimensions appear in 
the OPE comparing to the QCD sum rules in vacuum. Thermal version of QCD sum rules 
has been successfully used to study the thermal properties of light [§-11], heavy-light H-14 
and heavy-heavy [15l-ll8| mesons as a trusty and well-established approach. 

In the present work, we calculate the mases and decay constants of the heavy vector 
quarkonia J/ip (cc) and T (bb) in the framework of the thermal QCD sum rules. Using the 
thermal quark propagator, we calculate the expression for the spectral density in one loop 
approach. Taking into account also the two loop perturbative contributions in a s order [l|, y[ 
as well as new nonperturbative contributions arising in thermal QCD in addition to vacuum 
version, we acquire thermal QCD sum rules for the masses and decay constants. Using 
the results of the energy density for the interval T = (0 — 170) MeV obtained via Chiral 
perturbation theory [l9| as well as the values of the energy density and gluon condensates 



obtained in the region T = (100 - 170) MeV via lattice QCD [20|, |2j, we present the 
sensitivity of the masses and decay constants of the J / ip and T heavy vector mesons on the 
temperature. In our calculations, we also use the temperature dependent two loop expression 
for the strong coupling constant obtained using the perturbation theory and improved by the 
lattice results [l6|, 22[. We see that the values of the decay constants decrease considerably 



near to the critical or deconfinement temperature comparing to their values in vacuum. This 
can be considered as a sign of the QGP phase transition. 

The rest of the paper is organized as follows. In the next two sections, we derive thermal 
QCD sum rules for the considered observables. The last section is devoted to the numer- 
ical analysis of the observables and present their temperature dependency as well as our 
discussion. 
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II. OPE OF THERMAL CORRELATION FUNCTION FOR HEAVY-HEAVY 

VECTOR MESONS 

To obtain the thermal QCD sum rules for physical quantities, we need to calculate the 
convenient thermal correlation function in two different ways: in terms of QCD degrees 
of freedom and in terms of hadronic parameters. In QCD side, the correlation function is 
evaluated via OPE which helps us expand the time ordering product of currents in terms 
of operators with different dimensions. In the present section, we obtain the OPE for the 
considered quantities. We begin by considering the following two point thermal correlation 
function: 

IV(g,T) = ijd'x <j**Tr(p r(j M (x)jt(0))), (1) 

where =: Q(x) , -f fl Q(x) : with Q = b or c is the vector current, T indicates the time 

ordered product and p = /Tre~^ H is the thermal density matrix of QCD at temperature 
T = 1//3. As we previously mentioned, the Lorentz invariance breaks down via the choice 
of reference frame at which the matter is at rest. However, using the four velocity vector 
of the matter, we can define Lorentz invariant quantities such as u = u ■ q and q 2 = u 2 — q 2 . 
By the help of these quantities, the aforesaid thermal correlation function can be expressed 
in terms of two independent tensors P^ u and at finite temperature i.e., 

n^(g,T) =g^n i (g 2 ,^) + P^n 4 (g 2 ,o;), (2) 

where 

p I Ml - it- ~ 

Q [11/ (^) 

and Up = Up — u q^j q 2 . The functions IT; and Il t are the following Lorentz invariant functions: 

rb(V,u;) = l^n^, (4) 



U t (q 2 ,co) = -U g rU^ + ^U liV u v ). (5) 



t 
t 

It can be shown that in the limit |q| — > 0, the lb function can be expressed as lit 



— ^g^H^ [10] and one can easily find the H-tylo, |q| = o) = q% rLYg , |q| = o) relation 

between two II; and IT t functions. In real time thermal field theory, the function Ui(q 2 ,u) 
or U t (q 2 ,oj) can be written in 2 x 2 matrix form and elements of this matrix depend on 
only one analytic function [23]. Therefore, calculation of the 11-component of this matrix is 
sufficient to determine completely the dynamics of the corresponding two-point function. It 
can also be shown that in the fixed value of |q|, the spectral representation of the thermal 
correlation function can be written as [9|]: 

TT ( 2 rp\ r , ,2 Pl,t(q'o 2 ,T) 

TkM,T)= dq ,2 ; n2 L , (6) 
v ' J o q +Qo 



4 



where Ql 



-Qq, and 



pi, t (q%,T^ = -ImTi ht (ql,T} tanh^. 



7T ^ ' 2 

The thermal correlation function of Eq. ([1]) can be written in momentum space as: 

d 4 k 



n^(g,T) 



(27T) 



-Tr 



(7) 



(8) 



where, we consider the 11-component of the S(k) (thermal quark propagator) which is 
expressed as a sum of its vacuum expression and a term depending on the Fermi distribution 
function ]2M 



S{k) = ( 7 % + m)(- l - r — + 2mn(\h\)5(k 2 - m 2 )) 

\h z — m + IE 



(9) 



where n(x) = [exp(/3x) + is the Fermi distribution function. Now, we insert the prop- 
agator of Eq. (|9]) in Eq. ([8]) and consider IIi(g, T) = g^ILj^g, T) function. Carrying out 
the integral over k Q , we obtain the imaginary part of the IIi(g, T) in the following form: 



ImIL 1 (q,T) = L{q ) + L(-q ), 



(10) 



where 
L(q ) = N t 



<ik cu 2 — k + k ■ q — uiqo — 2m 2 



An 2 

x ([(1 - m)(l - n 2 ) + nin 2 ]5(qo -u x - u 2 ) - [(1 - ni)n 2 + (1 - n 2 )ni]8(q -u t + w 2 )), 

(11) 



and n\ = n(iOx), n 2 = n(i0 2 ), coi = \J~k 2 + m 2 and co 2 = \J(k — q) 2 + m 2 . The terms without 
the Fermi distribution functions show the vacuum contributions but those including the 
Fermi distribution functions depict medium contributions. The delta-functions in the differ- 
ent terms of Eq. ffTTj) control the regions of non- vanishing imaginary parts of Hi(q, T), which 
define the position of branch cuts [5(. After straightforward calculations, the annihilation 
and scattering parts of pi(q 2 ,T^J = ^ImTli(q 2 , T^j tanh ^ at nonzero momentum can be 



written as: 
Pl,a 



-3g 2 



57T 



(3 - v 2 ) 



dx n + (x) 



for 



4m 2 + q 2 < q 2 < 00, 



(12) 



Ply 



3q 2 

167T : 



■(3 



v 



dx 



n~(x) — n+(x) 



for q 2 < q 2 , 



n 



(13) 



2 (|q|a; — go) • From 



where z^(g 2 ) = \Jl — 4m 2 /g 2 , n+(x) = n |(go + |q|^) and n-(x) 
the similar manner, one can calculate also the function H 2 (q, T) = uTI^g, T)u v . Using the 
obtained results in Eqs. (jlj) and (JSJ), the annihilation and scattering parts of p t at nonzero 
momentum is obtained as: 

3g 2 

Pt,a ~ 327r2 



dx(2 - u 2 + x 2 )[l - 2n+{x)}. 

J —v 



(14) 
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Pt,s = -^tt^ / dx ( 2 -v 2 + x 2 )[n_(x) - n+(x)]. 
62,11 Jv 



(15) 



The annihilation part of pi, i.e., p^ a and its scattering part p^ s also at nonzero momentum 
can be found from Eqs. ( TH| and (1151) replacing the coefficient (2 — v 2 + x 2 ) by 2(1 — x 2 ). 

In our calculations, we also take into account the perturbative two-loop order a* correction 
to the spectral density. This correction at zero temperature can be written as [H, 0] : 



p Qs (s) = a s -^v(s)(3 - is 2 (s) 



71 



2v(s) 



(3 



i/( s 



An' 



(16) 



where, we replace the strong coupling a s in Eq. (fT6"j) with its temperature dependent lattice 
improved expression a(T) = 2.095(82)^ HQ, here 



8tt 2 VA- 



21n(- ) 



(17) 



where A^ = T c /1.14(4) and T c = 0.160GeV. 

Now, we proceed to calculate the nonperturbative part in QCD side. For this aim, we use 
the nonperturbative part of the quark propagator in an external gluon field, A^(x) in the 
Fock-Schwinger gauge, x^A^x) = 0. Taking into account one and two gluon lines attached 
to the quark line, the massive quark propagator can be written in momentum space as jsj]: 



gaa'nonpert^ = * (^aa' gc^ty 



where, 



(k 2 — m 2 y 



cr K \{ft + m) + (ft + m)a K \ 



ft + m 
(k 2 — m 2 ) 5 



[faf3[iu + fa^v + fanup){ft + Til), 



faPtiv = la(ft + m)jp(ft + mjj^ift + m)j u . 



(18) 



(19) 



To go on, we also need to know the expectation value (TrGa^G^) . The Lorentz covariance 
at finite temperature permits us to write the general structure of this expectation value in 
the following manner: 



{Tr c G a pG^ p ) — —(g afJ ,9i3u — 9au9/3ii)(GxaG a a ) 



1 

+ 6 



g af i9fi v - 9 a u9p^ - 2{uaU^gpv - u a u u gp^ - u p u^g au + UpU v g atl ) (u x Q 9 Xa u a ), 

(20) 



where, as we also previously mentioned is the four-velocity of the heat bath and it is 
introduced to restore Lorentz invariance formally in the thermal field theory. In the rest 
frame of the heat bath «'* = (!, 0, 0, 0) and u 2 = 1. Furthermore, Q g Xa is the traceless gluonic 
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part of the stress-tensor of the QCD. Up to terms required for our calculations, the non 
perturbative part of massive quark propagator at finite temperature is obtained as: 



1 1 

Qaa 1 nonpert ( i\ „/iC\ao'/ic 
' ~V I{ ' ] G "*(]fc2- m 2)2 

J2 raa' 



cr KA (f + m) + (ft + m)a K \ 



i g 2 5 aa r m[k 2 +m jf) 



1 



3 {k 2 -m 2 )^ 4 V afi ' 3(k 2 -m 2 



m{k 2 - m 2 )(k 2 - A{k ■ u) 2 ) 
+ (m 2 -k 2 )(-m 2 + A(k-u) 2 ) ft + A(k ■ u)(m 2 - k 2 ) 2 i]{u a Q 9 ap u p )}. (21) 



Using the above expression and after straightforward but lengthy calculations, the nonper- 
turbative part in QCD side is obtained as: 



n 



nonpert 



(a s G 2 ) 



72n 



m 2 + q 2 (— 1 + x)x 



6g b (-l + x) V + 6m 2 gV(-l + x) 2 (l - Qx + Qx 2 



+ m 6 (5 - 32a; + 42x 2 - 20x 3 + 10x 4 ) + m A q 2 x( K - 14 + 95x - 140x 2 + 65 3 + 6x 4 - 2x 



547T 



m 2 + q 2 (— 1 + x)x 



x(-l + x)(4gV(l - 3x + 2x 2 ) 2 + m 4 (12 - 35x + 21x 2 + 28x 3 



- 14a; 4 ) + mVa;(-13 + 55x - 82a; 2 + 36a; 3 + 6a; 4 - 2a; 5 )) (q 2 - 4{q ■ u) 2 )] }, (22) 
where, (G 2 ) = (G c a pG ca ?). 



III. PHENOMENOLOGICAL PART AND THERMAL SUM RULES 

Now, we turn our attention to calculate the physical or phenomenological side of the 
correlation function. For this aim, we insert a complete set of physical intermediate state to 
Eq. ([I]) and perform integral over x. Isolating the ground state, we get 

E (W»M*) + 

a m v ~ Q 

where the hadronic states {|U(g, A))} form a complete set and indicate the contributions 

of excited vector mesons and continuum states. 

In order to obtain thermal sum rules, now we equate the spectral representation and 
results of operator product expansion for amplitudes IU(g 2 ,u;) or Ib(g 2 ,(X>) at sufficiently 
high Qq. When performing numerical results, we should exchange our reference to one at 
which the particle is at rest, i.e., we shall set |q| — > 0. In this limit since the functions II; and 
IT t are related to each other, it is enough to use one of them to acquire thermal sum rules. 
Here, we use the function Yl t . When we use the standard spectral representation, if the 
spectral density at s — > oo does not approach to zero, in this case the correlation function 
is expressed in terms of a diverge integral. In such a case, to overcome this problem, we 
subtract first few terms of its Taylor expansion at q 2 = from rb(g 2 ,u;), 

H<* 1,1) - n,,o) + Q^Ql + f f (24) 
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Equating the OPE and hadronic representations of the correlation function and applying 
quark-hadron duality, our sum-rule takes the form: 

mi _ Q 4 r hM±^A ds + /' q|2 _i±±_ ds + nr— , (25) 



(m 2 / + Qg) ^ u s 2 (s + Qg) Jo s + Q 2 

where, for simplicity, the total decay width of meson has been neglected. The decay constant 
fv is defined by the matrix element of the current between the vacuum and the vector- 
meson state, i.e., 

(0\J,\V(q,\)) = fvm v e^. (26) 

In derivation of Eq. (I25p we have also used summation over polarization states, J2\ £ p * £ 1 = 
~{9^ v ~ lulu/my)- The Borel transformation removes subtraction terms in the dispersion 
relation and also exponentially suppresses the contributions coming from the excited reso- 
nances and continuum states heavier than considered vector ground states. Applying Borel 
transformation with respect to Q$ to both sides of Eq. f[2"51 ) , we obtain 

f 2 m 2 v exp ( - ^0 = J 2 rfs [pt, a ( s ) + Pa a (s)}e~^ + ds p t , s (s)e~^ + Bn™ n P ert . 

(27) 

As we also previously mentioned, when doing numerical analysis, we will set |q| — > repre- 
senting the rest frame of the particle. In this case, the scattering cut shrinks to a point and 
the spectral density becomes a singular function. Hence, the second term in the right side 
of Eq. ( 12T1) must be detailed analyzed. Similar analysis has been also performed in (5], @|. 
Detailed analysis shows that 

Hm^V^^--^)^. (28) 

In Eq. ( 12 7p . B]j™ on P ert shows the nonperturbative part of QCD side in Borel transformed 
scheme, which is given by: 

^nonpert = f 1 j [ - ^ -1 ( (a s G 2 ) [l2 M & X* (-1 + xf 

1 Jo 144 Ti M 6 x (—1 + x) P Lm 2 x (-1 + x)J^ 7 L v ; 

- m 6 (1 - 2x) 2 (-l - x + x 2 ) - 12 m 2 M 4 x 2 (-1 + x) 2 (1 - 3x + 3x 2 ) + m A M 2 x (-2 + 19x 

- 32x 2 + llx 3 + 6x 4 - 2x 5 )] + 4 a s (G 9 ) \ — 8 M 6 x 3 (1 — 2x) 2 (-l + x) 3 + m 6 (1 - 2x) 2 (-l - x 
+ x 2 )-2 m 2 M 4 x 2 (-1 + a;) 2 (-l - Qx + 8x 2 - Ax 3 + 2x 4 ) + m 4 M 2 x (-2 + 3x - 12x 2 + 31x 3 
-30x 4 + 10x 5 )]}, (29) 

where, Q 9 = 6q . 



IV. NUMERICAL ANALYSIS 

In this section, we discuss the sensitivity of the masses and leptonic decay constants of 
the J/ip and T vector mesons to temperature and obtain the numerical results for these 
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quantities in vacuum. Taking into account the Eqs. (128"]) and f|29|) and applying derivative 
with respect to 1/M 2 to both sides of the Eq. (|27|) and dividing by themselves, we obtain 



m 2 v (T) 



f^P ds s [ptA s ) + p<*M\ ex P ( - jfr) 


+ nr nper '(M 2 ,T) 


C ( 2 T) ds\pt, a (s) + p as (s)} exp ( - j^) 


| ^-j-j-nonperi 



where 



and 



UT nper \M 2 ,T) = M 4 



Pt,a( S ) 



; su(s)(3-u 2 (s)) 



\2T , 



(30) 



(31) 



(32) 



As we did also i n [l8j| . we use the gluonic part of the energy density both obtained from 
lattice QCD [20|, l2l| and chiral perturbation theory [l9j . In the rest frame of the heat bath, 
the results of some quantities obtained using lattice QCD in [20] are well fitted by the help 
of the following parametrization for the thermal average of total energy density (0): 



(0) = 2(0 9 ) = 6 x 10- 6 ea:p[80(T - 0.1)](GeV 4 



(33) 



where temperature T is measured in units of GeV and this parametrization is valid only 
in the region 0.1 GeV < T < 0.17 GeV. Here, we should stress that the total energy 
density has been calculated for T > in chiral perturbation theory, while this quantity has 
only been obtained for T > 100 MeV in lattice QCD (see [13, 21] for more details). In low 
temperature chiral perturbation limit, the thermal average of the energy density is expressed 
as 0: 



(©) = (m + s P , 



(34) 



where (0p is trace of the total energy momentum tensor and p is pressure. These quantities 
are given by: 



7T 



rpS 



V 



270 F£ 

'm n T 



T 



3T(— ) 2 (1 + + 



15 T 105 T 2 



128 ml 



exp(^ 



T 



(35) 



where A p = 0.275 GeV, F w = 0.093 GeV and = 0.14 GeV. 

The next step is to present the temperature dependent continuum threshold and gluon 
condensate. In the present work, we use the sq(T) [l3[ and (G 2 ) 20, 21] as: 



where T* = 1.1 T r 



so(T) 
0.176 GeV. 

(G 2 ) 



so 



(t*) 



(0\G 2 \0) 



exp 



12(£ 



1.05 



(36) 



(37) 
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FIG. 1. The dependence of the mass of J/ip meson in vacuum on the Borel parameter M 2 



10 



0.6 



0.4- 



10 



15 



15 

M 2 (GeV 2 ) 



20 



25 



s 


=11 GeV 2 


..... 


=12 GeV 2 


— s 


=13 GeV 2 




0.6 



0.4 



20 



25 



FIG. 2. The dependence of the leptonic decay constant of J ftp meson in vacuum on the Borel 
parameter M 2 . 

In further analysis, we use the values, m c = (1.3 ± 0.05) GeV, m b = (4.7 ± 0.1) GeV 
and (0 | ±a s G 2 | 0) = (0.012 ± 0.004) GeV 4 for quarks masses and gluon condensate at 
zero temperature. The sum rules for the masses and decay constants also include two more 
auxiliary parameters: continuum threshold So and Borel mass parameter M 2 . These are not 
physical quantities, hence the physical observables should be approximately insensitive to 
these parameters. Therefore, we look for working regions of these parameters such that the 
dependences of the masses and decay constants on these parameters are weak. The contin- 
uum threshold, s is not completely arbitrary, but it is related to the energy of the first exited 
state with the same quantum numbers as the interpolating currents. Our numerical analysis 
show that in the intervals so = (11 — 13) GeV 2 and so = (98 — 102) GeV 2 , respectively for 
the J/ip and T channels, the results weakly depend on this parameter. The working region 
for the Borel mass parameter, M 2 is determined demanding that both the contributions of 
the higher states and continuum are sufficiently suppressed and the contributions coming 
from the higher dimensional operators are small. As a result, the working region for the 
Borel parameter is found to be 8 GeV 2 < M 2 < 25 GeV 2 and 12 GeV 2 < M 2 < 35 GeV 2 
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FIG. 3. The dependence of the mass of T meson in vacuum on the Borel parameter M z 
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FIG. 4. The dependence of the leptonic decay constant of T meson in vacuum on the Borel 
parameter M 2 . 



in J/if> and T channels, respectively. 

Using the working regions for auxiliary parameters as well as other input parameters, we 
depict the dependence of the masses and leptonic decay constants of the heavy J /if) and T 
vector quarkonia in Figs. (JTJSI) at T = (vacuum). By a quick glance in these figures, we see 
that the masses and decay constants represent good stability with respect to the variation 
of the Borel parameter in its working region. Also, we see a weak dependence of the results 
on the continuum threshold, sq. From these figures, we deduce the numerical values of these 
parameters as shown in the Tables (jl]) and (JTTJ) . The uncertainties presented in these Tables 
are due to the uncertainties in calculation of the working regions for the auxiliary parameters 
as well as errors in the values of the other input parameters. In these Tables, we also present 
the existing numerical predictions from the other approaches such as lattice QCD, potential 
model and nonrelativistic quark model as well as existing experimental data. As far as 
the leptonic decay constants are concerned, our predictions for the central values are a 
bit bigger than the predictions of the other approaches and experiment, but when taking 
into account the uncertainties, our results become comparable especially with the potential 
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fj/^MeV) 


f r (MeV) 


Present Work 


481 ± 36 


746 ± 62 


Lattice [25, 26] 


399 ± 4 




Experimental [25, 261 


409 ± 15 


708 ±8 


Potential Model [25] 


400 ± 45 


685 ± 30 


Nonrelativistic Quark Model [26] 


423 


716 



TABLE I. Values of the leptonic decay constants of the heavy-heavy J/ip and T vector mesons in 
vacuum. 





mj/il> {GeV) 


m r (GeV) 


Present Work 


3.05 ± 0.08 


9.68 ±0.25 


Experimental [27] 


3.096916 ±0.000011 


9.46030 ± 0.00026 



TABLE II. Values of the masses of the heavy-heavy J/ip and T vector mesons in vacuum. 

0.00 0.05 0.10 0.15 0.20 

3.4- 

3.2- 



""i 3.0 -, 
E 

2.8- 
2.6- 



0.00 0.05 0.10 0.15 0.20 

T(GeV) 



FIG. 5. The dependence of the mass of J/ip vector meson in GeV on temperature at M 2 = 10 GeV 2 . 

and nonrelativistic quark models predictions as well as experimental data. However, our 
predictions on the masses are in good consistency with the experimental values. 

Our final task is to discuss the temperature dependence of the leptonic decay constant 
and masses of the considered particles. For this aim, we plot these quantities in terms of 
temperature in figures ([ME]) using the total energy density from both chiral perturbation 
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FIG. 6. The dependence of the leptonic decay constant of J/ip vector meson in GeV on temperature 
at M 2 = 10 GeV 2 . 
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FIG. 7. The dependence of the mass of T vector meson in GeV on temperature at M 2 = 20 GeV 2 . 

theory and lattice QCD (valid only for T > 100 MeV) and at different fixed values of the 
So but a fixed value of the Borel mass parameter. From these figures, we observe that the 
masses and decay constants remain insensitive to the variation of the temperature up to 
T = 100 MeV, however after this point, they start to diminish increasing the temperature. 
At deconfinement or critical temperature, the decay constants approach roughly to 45% of 
their values at zero temperature, while the masses are decreased about 12%, and 2.5% for 
J/ip and T states, respectively. Considerable decreasing in the values of the decay constants 
near the deconfinement temperature can be judged as a sign of the quark gluon plasma 
phase transition. 
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FIG. 8. The dependence of the leptonic decay constant of T vector meson in GeV on temperature 
at M 2 = 20 GeV 2 . 
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